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The structural correlates of functional dysconnectivity in autism spectrum disorders (ASD) 
have been seldom explored, despite the fact that altered functional connectivity is 
one of the most frequent neuropathological observations in the disorder. We analyzed 
cerebral morphometry and structural connectivity using multi-modal imaging for 11 
children/adolescents with ASD and 11 matched controls. We estimated regional cortical 
and white matter volumes, as well as vertex-wise measures of cortical thickness and local 
Gyrification Index (/Gl). Diffusion Tensor Images (DTI) were used to measure Fractional 
Anisotropy (FA) and tractography estimates of short- and long-range connectivity. We 
observed four clusters of /Gl reduction in patients with ASD, three were located in the 
right inferior frontal region extending to the inferior parietal lobe, and one was in the 
right medial parieto-occipital region. Reduced volume was found in the anterior corpus 
callosum, along with fewer inter-hemispheric frontal streamlines. Despite the spatial 
correspondence of decreased gyrification and reduced long connectivity, we did not 
observe any significant relationship between the two. However, a positive correlation 
between /Gl and local connectivity was present in all four clusters in patients with ASD. 
Reduced gyrification in the inferior fronto-parietal and posterior medial cortical regions 
lends support for early-disrupted cortical growth in both the mirror neuron system and 
midline structures responsible for social cognition. Early impaired neurodevelopment in 
these regions may represent an initial substrate for altered maturation in the cerebral 
networks that support complex social skills. We also demonstrate that gyrification 
changes are related to connectivity. This supports the idea that an imbalance between 
short- and long-range white matter tracts not only impairs the integration of information 
from multiple neural systems, but also alters the shape of the brain early on in autism. 
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INTRODUCTION 

Autism is a heterogeneous disorder characterized by a triad of 
symptoms including impairments in social interactions, delayed 
development of spoken language, and repetitive patterns of 
behavior. To satisfactorily account for the observed clinical 
heterogeneity in autism, the name "autism spectrum disor- 
der" (ASD) is commonly used to convey the associated clini- 
cal manifestations that vary in severity along a continuum of 
autistic traits. Most epidemiological records estimate the global 
prevalence of ASD at 1 in 160 individuals (Elsabbagh et al, 
2012), with some studies reporting rates as high as 1 in 88 
children (Centers for Disease Control and Prevention, 2008). 
Understanding the neurobiological bases of this pervasive devel- 
opmental disorder, which highly impacts the societal integra- 
tion and professional achievements of affected persons, is one 



of the main motivations driving the prolific research on the 
disorder. 

Structural and functional neuroimaging studies on ASD are 
particularly abundant. Early morphometric studies reported 
increased brain volume, but decreased total volume thereafter, 
during the first years of life in patients with ASD compared to 
healthy controls (Courchesne et al, 2001; Courchesne, 2004). 
This pattern of early overgrowth followed by degenerative change 
has been explained by a failure to refine the cerebral cir- 
cuitry through the adaptive, experience-driven processes nor- 
mally occurring during childhood (Courchesne et al, 2011). 
Indeed, there is a large amount of evidence for disrupted 
organization of cerebral networks in children, adolescents and 
adults with autism. Structural brain imaging studies reported 
increased white matter volume in regions corresponding to local 
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cortico- cortical connections, and only minor changes, or even 
sometimes decreased volumes, in regions corresponding to long- 
range or inter-hemispheric connections (reviewed in Minshew 
and Williams, 2007). Increased local connectivity and reduced 
distant connectivity was further corroborated by post-mortem 
examinations (Zikopoulos and Barbas, 2010) and by functional 
studies using EEG (Barttfeld et al, 2011). Finally, several fMRI 
studies measuring functional connectivity using resting- state 
paradigms observed decreased long-range functional connectivity 
in children or young adults with ASD (Kennedy and Courchesne, 
2008; Assaf et al, 2010; Weng et al, 2010; von dem Hagen et al, 
2013). Despite the plethora of evidence for disrupted structural 
and functional connectivity and a growing body of research 
demonstrating morphometric differences in the brains of patients 
with autism; there are surprisingly few integrated studies that 
show how differences in cerebral morphology and connectivity fit 
together. Adding to our knowledge of the relationships between 
different anatomical impairments should elucidate the mecha- 
nisms underlying brain alterations in autism. 

It is accepted that cortical folding reflects a person s prenatal 
development (Regis et al., 2005) (as well as events from the 
first months of post-natal life (Schaer et al., 2009; Haukvik 
et al, 2011)). It follows that measuring the shape of the cor- 
tex, using three-dimensional cortical reconstructions, provides us 
with insight into early brain development. Although the processes 
underlying the creation of specific sulcal patterns are poorly 
understood, existing theories point to the determinants of early 
cortical folding. Initial hypotheses proposed that gyrification 
results from mechanical forces intrinsic to the cortex (Richman 
et al., 1975; Welker et al., 1990). However, more recent theories 
consider cortical shape as a product of underlying patterns of 
connectivity, implicating alterations to both connectivity and 
cortical folding, both of which are highly relevant in autism. The 
tension-based model of convolutional development (Van Essen, 
1997) postulates that strongly interconnected cortical regions are 
pulled towards one another during embryological development, 
resulting in both compact and streamlined wiring of the nervous 
system. According to this model, disturbed gyrification in the 
adult human brain reflects abnormal patterns of white matter 
connectivity. Measuring gyrification abnormalities at any age may 
signal early adverse events and contribute to our understanding of 
both the timing and the nature of brain alterations in neurodevel- 
opmental disorders. 

Previous studies have noted alterations to cortical shape in 
autism (Levitt et al, 2003; Nordahl et al, 2007; Shokouhi et al, 
2012), and some have used the Gyrification Index (GI; Hardan 
et al, 2004; Casanova et al, 2009; Kates et al, 2009; Jou et al, 
2010; Meguid et al., 2010). Given that the cortex grows primarily 
through radial expansion (Rakic, 1988), the GI was specifically 
designed to identify early defects in cortical development. How- 
ever, all but one (Wallace et al., 2013) existing studies quanti- 
fying GI in patients with ASD used two-dimensional or global 
approaches, which do not allow for the identification of focal 
differences. By contrast, the local Gyrification Index (/GI; Schaer 
et al., 2008) has been shown to provide reliable estimates of GI 
with fine-grained resolution in many conditions (Zhang et al., 
2009; Juranek and Salman, 2010; Zhang et al, 2010; Palaniyappan 



and Liddle, 2012; Palaniyappan et al., 2011; Ronan et al., 2011; 
Thesen et al., 2011; Srivastava et al., 2012). 

In the present study, we sought to combine advanced multi- 
modal techniques in a small individually- matched group of ASD 
and healthy controls in order to simultaneously examine alter- 
ations in gyrification and structural connectivity. We first propose 
an exploratory analysis of the morphometry of gray and white 
matter structure, and of white matter connectivity. For that pur- 
pose, we use the Tl -weighted imaging to quantify the total cere- 
bral and cerebellar gray and white matter volumes, the regional 
cortical and white matter volumes, and to measure cortical thick- 
ness and /GI at thousands of points across the hemisphere. We also 
exploit Diffusion Tensor Imaging (DTI) to quantify voxel- wise 
alterations in white matter microstructure and use tractography 
to provide estimates of structural connectivity. In subsequent 
analyses, we aim at integrating the morphometric and the con- 
nectivity findings. In line with the Van Essen's tension-based 
morphogenesis hypothesis, we expect to observe that regions with 
altered gyrification in autism corresponds to areas with aberrant 
patterns of short- and long-range connectivity, as quantified using 
tractography. 

MATERIALS AND METHODS 

PARTICIPANTS 

Patients with ASD 

Eleven children and adolescents with ASD participated in the 
current study (8 males). The group had an average age of 12.9 =b 
2.7 years (range 9.3-17.4) and an average full-scale IQ (using the 
Wechsler WISC-III (Wechsler, 1991)) score of 79.4 zb 18.1 (range 
51-105). Participants were recruited with the help of local asso- 
ciations, therapeutic schools, and a local ASD diagnostic clinic. 
Once participants contacted us to express interest in the study, a 
detailed medical history, including details about their diagnosis 
were taken. Individuals with known genetic disorders, as well 
as malformations and birth defects, were excluded. An initial 
appointment was then set to reconfirm participants' diagnoses 
using the autism diagnostic interview- revised (ADI-R) interview 
with one or both parents. The group of patients with ASD had 
the following scores on the ADI-R (Le Couteur et al., 2007): social 
interaction: 15.7 =b 7.9, communication: 12.4 zb 6.6, restricted and 
repetitive behaviors: 5.7 =b 3.3. The ADI-R was followed by an 
autism diagnostic observation scale (ADOS; Lord et al, 2009), 
which was conducted by a research- reliable clinician from the 
institution s ASD diagnostic clinic. The parents of all participants 
also filled in the Social Communication Questionnaire (SCQ; 
Berument et al, 1999). 

Control participants 

The comparison group was comprised of 11 healthy controls, 
individually matched with each patient for gender and age. 
The control group had an average age of 12.7 =b 2.7 (range 
8.7-16.8). There was no difference in mean age between patients 
and controls {p = 0.897). The average full-scale IQ of the control 
group was 110.5 =b 13.3 (range 88-129). 

Written informed consent was received from all subjects and 
their parents in accordance with protocols approved by the local 
ethics committee. 
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IMAGE PROCESSING 

Cerebral magnetic resonance images were acquired using a 
Siemens Trio 3T scanner at the Geneva Center for Biomedical 
Imaging (CIBM). A set of Tl -weighted 3D volumetric images was 
acquired as a series of 192 contiguous coronal slices, with a voxel 
size of 0.86 x 0.86 x 1.1 mm (repetition time (TR) = 1200 ms, 
echo time (TE) = 3 ms, flip angle = 8°). DTI were acquired on 
the same day as a series of 64 axial slices with 30 directions, with 
a voxel size of 2 x 2 x 2 mm (bO = 1000 ms, TR = 8300 ms, 
TE = 82 ms, flip angle = 90°). 

Cortical reconstruction 

The Tl -weighted images were used to create cortical reconstruc- 
tion and volumetric segmentation using the FreeSurfer package 
(Martinos Center for Biomedical Imaging, Massachusetts Gen- 
eral Hospital, Boston^). Briefly, the processing was comprised 
of removing non-brain tissue, executing automatic segmentation 
of the subcortical gray matter structures, and extracting cortical 
surface (Dale et al., 1999; Fischl, 2012). Both intensity and con- 
tinuity data from the entire three-dimensional MR volume are 
used in the segmentation procedures, thus producing accurate 
representations of cortical thickness and volumes. These pro- 
cedures have been validated against histological studies (Rosas 
et al., 2002) and have been shown to be reliable across scanner 
models and field strengths (Han et al, 2006). At the end of 
the reconstruction process, the following volumes were available: 
total cerebral gray and white matter volumes, cerebellar gray and 
white matter volumes, corpus callosum volume, and the volumes 
of subcortical structures including thalamus, putamen, pallidum, 
caudate nucleus, as well as amygdala and hippocampus. 

Regional cortical volumes 

Subsequent to cortical reconstruction, the cortex was also sub- 
divided into units based on gyral and sulcal structures (Desikan 
et al, 2006). This parcellation method has been shown to be 
both valid and reliable, with high intra- class correlation coeffi- 
cients between the manual and automated procedures for both 
cortical volume estimates and region boundaries. The parcellation 
produces 34 cortical regions subdivided into 1 1 frontal regions, 
9 temporal regions, 5 parietal regions, 4 occipital regions, 4 parts 
of the cingulate cortex, and one label for the insula. 

Regional white matter volumes 

The parcellation of the cortical gray matter was subsequently used 
to subdivide the underlying white matter as described in Salat 
et al. (2009), a Voronoi diagram was created in the white matter 
voxels based on the distance to the nearest parcellation label, using 
a distance constraint of 5 mm. As a result of this process, regional 
white matter volumes were available for each of the 34 regions 
corresponding to the aforementioned gyral labeling. 

The corpus callosum was also identified and subdivided into 
5 portions along its anteroposterior axis, according to procedures 
detailed in Rosas et al. (2010). The volume of the corpus callosum 
was measured for each of the 5 portions (anterior, mid- anterior. 



^http://surfer,nmr.mgh,harvard.edu 



center, mid-posterior and posterior) on a 5 mm lateral extent 
centered on the mid-sagittal place. 

Cortical thickness and cortical gyrification 

Cortical thickness was measured in the native space of the images, 
as the shortest distance between the white (gray- white boundary) 
and the pial (gray-CSF interface) surfaces. As a result, cortical 
thickness values with submillimeter accuracy were available at 
more than 150,000 different points over each hemisphere res- 
olution (Fischl and Dale, 2000). Finally, based on the outer 
cortical surface reconstruction (pial surface), IGl was measured 
at thousands of points across each hemisphere using previously 
validated algorithms (Schaer et al., 2008). IGl is a surface-based 
measurement of the degree of cortical folding that iteratively 
quantifies the amount of cortex buried within the sulcal folds in 
the surrounding circular region. 

Inter- subject comparison of the cortical thickness and gyri- 
fication values is achieved through spherical registration of the 
surfaces that minimizes metric distortion and allows for a highly 
reliable point-to-point comparison of cortical thickness between 
groups (Fischl et al., 1999). 

Tract-based spatial statistics of the white matter structure 

The DTI images were used for voxelwise statistical analysis of the 
Fractional Anisotropy (FA) using Tract- Based Spatial Statistics 
(TBSS; Smith et al, 2006), which is part of FSL software.^ First, FA 
images were created by fitting a tensor model to the raw diffusion 
data using algorithms embedded in the FDT toolbox, followed by 
skull stripping. As described in the original protocol (Smith et al., 
2006, 2007), subjects' FA data were then aligned into a common 
space using nonlinear registration. Next, the mean FA image was 
created and thinned to create a mean FA skeleton that represents 
the centers of all tracts common to the group. Each subject's 
aligned FA data were then projected onto this skeleton and the 
resulting data were fed into voxelwise cross-subject statistics. 

Tractography analyses 

To relate the cortical anatomy with the underlying architecture of 
white matter fibers, we used tools embedded in the Connectome 
Mapping Toolkit (Daducci et al., 2012).^ Briefly, registration 
between the Tl -weighted and DTI images was completed using 
the bbregister function of FreeSurfer. The DTI images were pro- 
cessed with Diffusion Toolkit software"* using the deterministic 
streamline algorithm (Mori et al, 1999) to obtain tractographic 
reconstruction of white matter bundles. 

In the present study, we used the number of streamlines to 
quantify two different aspects of connectivity. First, we measured 
the amount of fibers connecting the homologous lobe between 
each hemisphere. The inter-hemispheric fibers obtained by this 
method represent a simple way to define long-range connectivity 
without having to define an arbitrary length threshold. To obtain 
the inter-hemispheric frontal fibers, we retained all streamlines 
connecting cortical regions comprised in the frontal lobe, as 



^ http ://fsLfmrib. ox. ac.uk/fsl/fslwiki/ 

^ http ://www. connectomics.org/connectomemapper/ 

^ http ://trackvis. org/dtk/ 
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defined in the Desikan parcellation (Desikan et al., 2006). To 
select the streamlines corresponding to inter-hemispheric parietal 
streamlines, we repeated the process with cortical regions corre- 
sponding to the parietal lobe. Finally, streamlines connecting the 
temporal and occipital cortical regions were considered together 
for this analysis, given the small amount of inter-hemispheric 
fibers connecting these two lobes. As a result, three different 
variables summarizing one aspect of long-range connectivity were 
studied: the number of inter-hemispheric frontal, parietal, and 
temporo- occipital streamlines. Second, we measured the connec- 
tivity within each lobe (i.e., the number of streamlines connecting 
one lobe to itself) as an estimate of short- connectivity that also 
doesn't require any arbitrary length threshold. 

STATISTICAL ANALYSES 
Volumetric analyses 

We used ANCOVA to compare cerebral, cerebellar and subcortical 
volumes between groups, including age and gender as covariates. 
To identify potential regional cortical alterations, we subsequently 
applied a MANCOVA on the 34 gyral regions in each hemisphere 
by entering diagnosis as the fixed factor, and both age and 
gender as covariates. Potential changes in regional white matter 
volumes were examined by doing a MANCOVA on the five corpus 
callosum regions, and another MANCOVA on the 34 regional 
white matter volume in each hemisphere. All the MANCOVA 
were performed with diagnosis as the fixed factor and both age 
and gender as covariates. 

Cortical thickness and gyrification analyses 

The comparisons of cortical thickness and gyrification over the 
whole brain used the fsaverage template included in the FreeSurfer 
distribution. Cortical thickness maps were smoothed using a 
fijll width at half maximum (FWHM) kernel of 10 mm. As IGl 
is already smooth (the degree of smoothness in our IGl data 
corresponds to a smoothing kernel of 10 mm), the data were 
not additionally smoothed prior to statistical analyses. Statistical 
analyses employed a General Linear Model (GLM) to estimate the 
effect of diagnosis, age and gender on thickness or gyrification at 
each cortical point. Cortical thickness or gyrification changes with 
age were fitted using a linear model. All results were corrected for 
multiple comparisons using the Monte Carlo simulation at the 
cluster level at the corrected significance threshold ofp < 0.05. 

Tract-Based Spatial Statistics of the white matter structure 

TBSS voxel-wise analyses were carried out across subjects for 
each point of the common skeleton. As our population was 
comprised of children and adolescents, the mean FA volume 
provided by FMRIB software library (FSL) ("FMRIB58") based 
on 58 adult brains was not optimal. Therefore, we chose the 
recommended alternative, using the "most typical" subject in our 
sample, to process the statistics. Local FA differences between 
patients and controls were tested for significance using a GLM. 
The skeleton-based approach has the advantage of reducing the 
number of statistical tests performed by reducing the num- 
ber of voxels being compared. Nevertheless, we performed a 
permutation-based approach to control for "Family- Wise Error" 
(FWE; Nichols and Holmes, 2002). The options we used in 
the statistic TBSS pipeline were the most recommended ones: 



the Threshold-Free Cluster Enhancement (TFCE) option and a 
number of permutations at 500 (see the TBSS user guide on 
http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/TBSS). A post-hoc t-test was 
ultimately conducted, comparing individual measurements of 
axial and radial diffusivity in clusters of significant between-group 
FA differences. 

Tractography analyses 

We used ANCOVA to compare the total number of streamlines 
with the total number of inter-hemispheric streamlines between 
patients with ASD and controls, while correcting for age and 
gender. We then conducted two MANCOVA to quantify potential 
differences between groups for the number of inter-hemispheric 
streamlines (long-range connectivity) and the number of stream- 
lines connecting each lobe with itself (short-range connecti- 
vity), with age, gender and the total number of streamlines as 
covariates. 

Correlations between gyrification and connectivity 

We conducted partial correlations between the average IGl in each 
cluster of between-group differences and measures of short- and 
long-range connectivity, while correcting for the effects of age, 
gender and total number of fibers. For long-range connectivity, 
we correlated IGl for each cluster with the inter-hemispheric fibers 
corresponding to the lobe where the largest part of the cluster 
was located. For the short-range connectivity, we correlated IGl 
with the number of intra-lobar fibers in the lobe where the largest 
part of the cluster was located. These partial correlations were 
conducted separately in ASD and control groups. Given that these 
correlations were based on our a priori hypothesis postulating a 
relationship between gyrification and underlying connectivity, we 
did not correct for multiple comparisons. 

Correlations with the clinical phenotype 

Finally, we explored how the neuroanatomical differences 
observed between the two groups may be related to the clinical 
outcome. For that purpose, we conducted partial correlations 
between neuroanatomical variables and the scores obtained at the 
ADI-R and in the SCQ, correcting for age and gender. In addition, 
variables measuring the number of streamlines between regions 
of interest were also corrected for the total number of streamlines. 
Correlations with clinical phenotype were not corrected for mul- 
tiple comparisons. 

RESULTS 

VOLUMETRIC ANALYSES 

We did not observe any significant differences between the 
cerebral, cerebellar and subcortical volumes of the groups (all 
p > 0.386). Upon further examination of the 34 cortical parcel 
volumes, we did not detect any significant patterns of change (left: 
Wilks Lambda: 0.042, p = 0.612; right: Wilks Lambda: 0.028, 
p = 0.517), nor did we detect pattern differences in the 34 
subcortical white matter regions (left: Wilks Lambda: 0.001, 
p = 0.087; right: Wilks Lambda: 0.037, p = 0.585). Despite the 
trend for a significant pattern of between-group differences in the 
MANCOVA corresponding to the left white matter subregions, 
none of the individual subregions revealed any significant differ- 
ence in the post-hoc analysis. 
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Selective reduction of the most anterior section of the corpus callosum in ASD 




Specific decrease in the number of interhemispheric frontal streamlines 
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FIGURE 1 I Convergent evidence for disrupted connectivity at the level of 
inter-hemispheric frontal connectivity in ASD. (A) Exannple of subdivision 
of the five sections of the corpus callosum displayed on a mid-sagittal slice. 
(B) Boxplots depicting volumetric measurements for the five sections in the 
two groups. The p-values are extracted from the MANCOVA (with correction 



for age and gender). (C) Example of tractographic reconstruction where the 
inter-hemispheric fibers are subdivided into three groups according to the 
cortical regions that they connect. (D) Boxplots comparing the number of 
streamlines connecting homologous lobes. The p-values are extracted from 
the MANCOVA (correction for age, gender and total number of streamlines). 



Results from the corpus callosum analysis are depicted in 
Figures lA, B. We observed significant group differences among 
the five sub-regions of the corpus callosum (Wilks Lambda: 0.354, 
p = 0.007, F = 5.11), with a selective reduction in the most 
anterior part of the corpus callosum in the ASD group compared 
to controls {p = 0.030). These results remained significant when 
covarying for total intracranial volume or total white matter 
volume instead of age and gender, as well as for total white matter 
volume, age and gender. 

VERTEX-WISE ANALYSES 

We did not observe any significant differences in cortical thick- 
ness related to diagnosis. However, vertex-wise comparisons of 
gyrification revealed four clusters of significant IGl reduction in 
patients with ASD compared to controls that remained significant 
after correcting for multiple comparisons. As shown in Figure 2, 



the clusters were all located in the right hemisphere, in the 
inferior parietal region, the lower part of the precentral gyrus, 
the inferior frontal gyrus, and the medial parieto-occipital region 
(cuneus/precuneus) . 

TBSS ANALYSES 

We found eight clusters of decreased FA in patients with ASD as 
compared to controls. The largest cluster of difference was located 
in the anterior part of the corpus callosum. The remaining seven 
clusters were located in the right hemisphere, no cluster of FA 
difference was seen in the left hemisphere. Figure 3 further details 
the distribution, location and size of these clusters. When compar- 
ing axial and radial diffusivity measurements in the clusters where 
FA significantly differed between patients with ASD and controls, 
we observed a significant between-group difference for axial dif- 
fusivity in only one cluster (cluster H, patients: 9.44e~^ =b 6.73e~^, 
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controls: 8.96e~^ =b 4.59e~^, p = 0.03), whereas significantly 
decreased radial diffusivity was observed in all eight clusters in 
patients with ASD compared to controls (allp < 0.003). 

TRACTOGRAPHIC ANALYSES 

We did not observe any difference in the total number of stream- 
lines reconstructed from the DTI images (ASD: 103021 zb 10883, 
controls: 107331 zb 12769, p = 0.341). 

We observed a significant reduction in the total number of 
inter-hemispheric fibers in patients with ASD (ASD: 4293 =b 
1834, controls: 6276 =b 1771; p = 0.033, F = 5.405). Further- 
more, and as demonstrated in Figures IC, D, we also observed 
a significant between-group difference in the regional pattern of 
inter-hemispheric fibers (MANCOVA covarying out the effect of 
age, gender and total number of fibers: Wilks Lambda: 0.478, 
p = 0.010, F = 5.47), showing a selective reduction in the number 
of inter-hemispheric frontal fibers (p = 0.002), with a selec- 
tive reduction in the number of inter-hemispheric frontal fibers 
(p = 0.002), but no significant differences in the inter- 
hemispheric parietal and temporo-occipital fibers. 

No difference in the pattern of short-range connectivity was 
observed between the group of patients with ASD and controls 
(Wilks Lambda: 0.600, p = 0.594). 

CORRELATIONS BETWEEN GYRIFICATION AND CONNECTIVITY 

Examining the relationship between gyrification and connectivity, 
we did not observe any significant relationships between the 
three clusters located in the frontal lobe and the number of 
inter-hemispheric frontal streamlines or in the occipital cluster 
and the number of inter-hemispheric occipital streamlines, in 



either diagnostic group. However, we found positive correlations 
between IGl and the variables measuring short-range connectiv- 
ity, in the ASD group only. As depicted in Figure 4, we observed 
a significant positive correlation between IGl in all three right 
frontal clusters and the number of streamlines connecting the 
right frontal lobe to itself {p = 0.043, p = 0.030 and p = 0.004 for 
clusters 1, 2 and 3 respectively as numbered on Figure 2). We also 
observed a positive correlation between IGl in the right occipital 
cluster (cluster 4 on Figure 2) and the number of streamlines 
connecting the right occipital lobe to itself {p = 0.010, R = 0.836). 
The spatial correspondence of the positive correlation between 
gyrification and short-range connectivity is further supported by 
the absence of a significant correlation between the three frontal 
clusters of IGl differences and the occipital connectivity, and the 
absence of a significant correlation between the occipital IGl and 
the frontal connectivity (as depicted in Figure 4 with dashed 
lines). 

CORRELATIONS WITH THE CLINICAL PHENOTYPE WITHIN ASD 

Results of the exploratory correlations between clinical scores at 
the ADI and SCQ and all variables that showed between-group 
differences are presented in Table 1. Briefly, IGl in the posterior 
cluster (number 4 in Figure 2) negatively correlated with the 
total score obtained on the ADI-R (R = -0.737, p = 0.024), 
and the reciprocal social interaction (R = —0.726, p = 0.027) 
and communication (R = —0.773, p = 0.015) domains from the 
ADI-R. Furthermore, the ADI-R scores in the domain of restric- 
tive and stereotyped patterns of behaviors negative correlated neg- 
atively with the number of inter-hemispheric fibers (R = —0.784, 
p = 0.021) and with the number of inter-hemispheric frontal 
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Table 1 | Correlations between clinical scores and neuroanatomical variables with between-group differences in participants with ASD. 



AD I: Total score 



ADI: Reciprocal 
social 
interactions 
score 



ADI: 
Communication 
score 



ADI: Restrictive 
and stereotyped 
patterns of 
behavior 



ADI: Early 
anomalies in 
development 



SCQ: Total score 



Gyrification 
IG\ in cluster 1 
/Gl in cluster 2 
/Gl in cluster 3 
/Gl in cluster 4 



-.005 (p = 0.989) 
-.136 (p = 0.728) 
-.037 (p = 0.925) 
-.737 (p = 0.024) 



Volumetric measurements 

Anterior section of .482 (p = 0.189) 

the corpus callosunn 



Tractograpliic measurements 



Nunnber of 
inter-hennispheric 
streannlines 
Nunnber of 
inter-hennispheric 
frontal streannlines 



-.345 (p = 0.402) 



-.390 (p = 0.339) 



-.216 (p = 0.577) 
-.386 (p = 0.305) 
-.292 (p = 0.446) 
.726 (p = 0.027) 

.537 (p = 0.136) 



.165 (p = 0.697) 



-.241 (p = 0.566) 



.030 (p = 0.940) 
-.020 (p = 0.959) 
.045 (p = 0.908) 
.773 (p = 0.015) 

.649 (p = 0.059) 



-.317 (p = 0.445) 
-.309 (p = 0.457) 



.497 (p = 0.173) 
.451 (p = 0.223) 
.494 (p = 0.177) 
-.043 (p = 0.912) 

-.325 (p = 0.394) 



-.784 (p = 0.021) 



-.779 (p = 0.023) 



.056 (p = 0.887) 
.041 (p = 0.917) 
.146 (p = 0.708) 
.205 (p = 0.597) 

-.419 (p = 0.262) 



.544 (p = 0.163) 



.392 (p = 0.337) 



.076 (p = 0.845) 
-.181 (p = 0.641) 
-.117 (p = 0.764) 
-.541 (p = 0.132) 

.315 (p = 0.409) 



-.682 (p = 0.063) 



-.681 (p = 0.063) 



This table provides R-values from partial correlations. Significance level is given in parentheses. Partial correlations accounted for an effect of age and gender on 
gyrification and volumetric measurements, as well as for an additional effect of total number of streamlines on tractographic measurements. Significant correlations 
at p < 0.05 (uncorrected) are highlighted in bold. 



fibers (R = —0.779, p = 0.023). None of the neuroanatomical 
variables correlated with the total score obtained on the SCQ. 

DISCUSSION 

In this study, we applied neuro imaging techniques using Tl- 
weighted and DTI images in the interest of quantifying morpho- 
metric and connectivity differences in a group of children and 
adolescents with ASD. We observed: (a) decreased gyrification in 
the right inferior frontal region extending into the inferior parietal 
region and in the medial parieto- occipital region of patients with 
ASD as compared to controls, the latter of which was related to the 
severity of social communications deficits in the group of ASD; 
(b) convergent evidence from three different analyses for altered 
long-range connectivity at the level of inter-hemispheric frontal 
fibers: volumetric reduction of the anterior corpus callosum, 
reduced FA in the anterior corpus callosum, and a decreased num- 
ber of virtual streamlines connecting homologous frontal lobes 
which further correlated with the severity of restrictive/repetitive 
behaviors; (c) further reduced FA in seven clusters of the right 
hemisphere of patients with ASD compared to controls; and (d) a 
positive correlation between /Gl in the clusters of between-group 
differences and short-range connectivity in the corresponding 
lobe. 

DECREASED GYRIFICATION 

We used a validated technique with exquisite resolution to mea- 
sure local cortical gyrification across the hemispheres (Schaer 
et al, 2008), and observed four clusters of reduced Gl in patients 
with ASD compared to controls, three of them located in the 
frontal lobe. This is in contrast with previous studies using Gl in 
children, adolescents or adults with ASD, which report either an 
increased Gl (Hardan et al., 2004; Jou et al., 2010), or an absence 



of significant difference (Casanova et al., 2009; Kates et al., 2009; 
Meguid et al, 2010). Lower intellectual abilities in our patient 
group may explain part of the divergence with previous results, 
given that both studies that reported increased Gl comprised 
participants with higher full-scale IQ scores (means: 105 =b 16 
for Hardan et al. (2004), 110 zb 15 for Jou et al. (2010) and 
113 zb 15 for Wallace et al. (2013)). However, we believe part 
of the difference to relate to the way Gl was calculated. Indeed, 
the two studies that reported higher Gl in the frontal lobe of 
subjects with ASD used manual delineation on one single frontal 
slice. Aside from the fact that manual tracing may be less reliable 
than automated delineation, measuring Gl on 2-D sections does 
not take into account the inherent 3-D nature of the cortical 
surface. 2-D measurement also can be biased by slice orientation 
(Zilles et al, 1997) and the presence of buried sulci (Magnotta 
et al., 1999), and it does not allow for precise localization of gyral 
anomalies in sublobar regions. Other studies that partly addressed 
these concerns did not report any significant differences in Gl in 
patients with ASD compared to controls (Casanova et al., 2009; 
Kates et al., 2009; Meguid et al., 2010). Casanova et al. used 
manual delineation in 40 randomly selected slices; Kates et al. 
applied an automated technique for measuring global and lobar 
Gl based on 2D sections; and Meguid et al. measured global Gl 
using three-dimensional cortical reconstructions. The technique 
that we use in the present study to measure IGl is automated, 
unbiased by slice orientation, and allows the quantification of 
gyrification differences at thousands of points over the recon- 
structed cortical surface. As a result, by using /Gl, we may have 
been able to detect gyrification differences in the frontal lobe of 
patients with ASD that had previously gone undetected. Studies 
using cortical reconstructions in autism corroborate this idea. 
They have detected focal changes to sulcal shape with alterations 
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Clusters of decreased FA in ASD compared to controls 



anteriori 



right 




FIGURE 3 I Differences in FA between patients with ASD and controls. 

Eight clusters of decreased FA in patients with ASD were observed at p < 
0.05 (corrected for multiple comparisons). Each cluster was attributed a letter 
according to decreasing size. The clusters were centered at the following 
MNI coordinates: A: x = 81 , y= 145, z = 90 (939 voxels); B: x = 52, / = 104, 



z = 104 (725 voxels); C: x = 66, / = 78, z = 101 (470 voxels); D: x = 57 K = 
106, z= 129 (186 voxels); E: x = 46, / = 131 , z = 112 (137 voxels); F: x = 82, 
)/ = 70, z = 93 (70 voxels); G: x = 43, / = 128, z = 90 (42 voxels); H: x = 60, 
/ = 65, z = 88 (23 voxels). We did not observe any clusters with increased FA 
in patients with ASD compared to controls. 



to the sylvian fissure and inferior fi-ontal sulcus (Levitt et al, 
2003), left fi-ontal operculum (Nordahl et al, 2007), and right 
intraparietal sulcus (Shokouhi et al., 2012). It is however worth 
noting that a recent study using the same technique as in the 
present study observed increased gyrification in different regions 
of the brain of 39 male adolescents with ASD as compared to 
41 controls, namely in bilateral occipital areas as well as in the 
left superior precuneus (Wallace et al., 2013). This discrepancy 
in the location and direction of gyrification changes using the 
same technique suggest either that different developmental mech- 
anisms take place in different regions of the brain of affected 
patients, or that that demographic characteristics (such as dif- 
ferences in age, gender, cognitive level, or symptom intensity) 
may have influenced the results. Indeed, it may be the case 
that the high clinical heterogeneity observed in patients affected 
with autism may be associated with different neurodevelopmental 
pathways. 

Decreased gyrification, as observed in the present study, is 
highly suggestive of reduced cortical expansion during early 
brain development, a process that might differentially affect spe- 
cific cortical regions. Neuropathological reports have pointed to 
abnormal cortical development in ASD, including a higher inci- 
dence of cortical dysgenesis, heterotopias and migration abnor- 
malities (Avino and Hutsler, 2010; Wegiel et al., 2010). Further 
detailed examination revealed that one cell type affected by migra- 
tion deficits in young children with ASD is von Economo neurons 
(Santos et al., 2011), which are spindle-shaped neurons thought 



to play a role in emotional function (Butti et al., 2013) that are 
located in the frontoinsular and cingulate cortices. The location of 
the von Economo neurons coincides with the location of cluster 
3 in the present study (see Figure 2). This anterior fronto -insular 
region is attracting increased attention in autism because of its 
key role in the "salience network" (Menon and Uddin, 2010). 
The anterior insula may have a critical role in processing infor- 
mation relevant to social functioning (Uddin and Menon, 2009) 
as a sort of "hub" that mediates interactions between cerebral 
networks that are processing information related to an external 
or internal stimulus. Functional neuroimaging studies tend to 
confirm the hypothesis of hypoactivation of the right anterior 
insula in autism, as pointed out by a meta-analysis based on 24 
functional neuroimaging studies examining social processes for a 
total of 276 patients with ASD and 291 controls (Di Martino et al., 
2009). 

Two other clusters, in the right inferior frontal gyrus (cluster 
2) and in a region extending from the right inferior part of the 
precentral gyrus to the inferior parietal region (cluster 1), are also 
located in regions that have received attention in ASD. Indeed, 
these regions are striking in their correspondence to the location 
of the fronto -parietal mirror neuron system, implicated in action 
imitation (Rizzolatti and Craighero, 2004). Decreased gyrification 
in the inferior fronto -parietal region thus supports altered devel- 
opment of the mirror system in ASD during in utero life or the 
first months after birth, pointing to a potential mechanism for 
early- disrupted abilities to imitate action of others. 
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FIGURE 4 I Correlations between gyrification and short-range 
connectivity within the group of patients with ASD. Partial 
correlations accounting for age, gender and total number of 
streannlines were conducted between /Gl in the cluster of 



between-group differences and the number of intra-lobar streamlines in 
the lobe where most of the cluster was located. Significant positive 
correlations were observed only in the corresponding lobe (dashed 
lines report non-significant relationships). 



The final cluster of reduced gyrification (cluster 4) is cen- 
tered in the occipital lobe, encompassing the cuneus and the 
pericalcarine sulcus, and further extending to the precuneus. 
Volumetric reductions have been consistently reported in the pre- 
cuneus in structural neuroimaging studies of ASD (Cauda et al, 
2011). The precuneus is a key part of the default mode network 
(DMN), which is thought to be concerned with self- referential 
and introspective activity, including the ability to understand 
others' intentions (Fair et al., 2008). Resting state paradigms have 
received increased recent interest in autism given the crucial role 
of the DMN in some aspects of social cognition. It is currently 
unclear to what extent the fronto-parietal mirror system (where 
we observed decreased /Gl in clusters 1 and 2) interacts with 
other regions of the social brain, including regions of the DMN. 
In an attempt to integrate these different views of the social 
brain, (Uddin et al, 2007) postulated that the cortical midline 
structures of the DMN and the fronto-parietal mirror-neuron 
system may represent two interwoven parts of self-related pro- 
cessing and social cognition: the mirror neurons encode physical 
aspects of social understanding (motor simulation and imitation 
of behaviors) and the midline DMN structures are associated 
with sophisticated processing of social interactions. Accordingly, 
reduced gyrification in the mirror-neuron system may impair 
physical aspects of the self- other relationship, consequently alter- 
ing the developmental cascade of the DMN responsible for more 



sophisticated social skills, such as empathy and theory of mind. 
The fact that cluster 4 extends into the precuneus may also point 
to an early defect, on top of which altered cortical maturation sub- 
sequently occurs. Indeed, the currently observed inverse relation- 
ship between gyrification and the level of autistic symptoms in 
the domains of reciprocal social interactions and communication 
points to the idea that early cortical development may determine 
the subsequent development of these more sophisticated social 
skills encoded in the precuneus. 

ABSENCE OF CORTICAL THICKNESS OR VOLUME DIFFERENCE 

We found reduced gyrification in the absence of differences in 
cortical thickness or volume. This absence contrasts with numer- 
ous studies that have reported altered cortical volume or thickness 
differences in ASD. Most studies have reported increased cortical 
volume or thickness in children with ASD (Hardan et al., 2006; 
Mak-Fan et al, 2012), whereas studies in adults have yielded 
more diverse results. Some studies in adults with ASD have 
reported mostly decreased cortical thickness (Hadjikhani et al., 
2006; Jiao et al., 2010; Wallace et al., 2010), some have shown 
a co-occurrence of thickening together with thinning (Ecker 
et al., 2010, 2013), and at least one has shown mostly thickening 
(Dziobek et al., 2010). Volumetric studies have more consis- 
tently reported increased volume in children and reduced volume 
in adults, supporting the hypothesis of early brain overgrowth 
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followed by neurodegenerative changes (Courchesne, 2004). 
Indeed, using a cross-sectional design with patients aged 1 to 
50, the largest study published to date provides evidence for an 
aberrant trajectory of cortical volume changes with age, with a 
pattern of early overgrowth during the first years of life, followed 
by decreased volume around 7 or 8 years old (Courchesne et al., 
201 1). Longitudinal studies confirm this pattern of abnormal cor- 
tical development in toddlers with ASD (Schumann et al, 2010) 
and of higher rates of cortical loss with age (Hardan et al., 2009). 
It should be noted, however, that studies recording abnormal 
trajectories of cortical features require large sample sizes, a broad 
age range at inclusion, and preferably, a longitudinal design. By 
contrast, our small study sample does not have the power to detect 
subtle cortical thickness differences that may be further diluted by 
complex maturational changes. 

CONVERGENT EVIDENCE FOR ALTERED LONG-RANGE CONNECTIVITY, 
MOST PROMINENT IN THE FRONTAL REGION 

As the largest white matter bundle of the brain, the corpus 
callosum represents the most essential component of connectivity, 
and more specifically long-range connectivity. Several fMRI and 
EEG studies have reported decreased long-range connectivity 
(reviewed in Belmonte et al., 2004). Patients with ASD were also 
shown to perform poorly on tests of inter-hemispheric transfer 
for auditory, visual and motor tasks (Nyden et al, 2004). More 
generally, decreased abilities in associative (Nikolaenko, 2001a) 
and metaphoric thinking (Nikolaenko, 2001b) were thought to 
depend on decreased inter-hemispheric information transfer. 
Here, we observed decreased volume of the most anterior part 
of the corpus callosum, reduced inter-hemispheric frontal con- 
nections, and decreased FA in the anterior corpus callosum, pro- 
viding strong multimodal evidence for altered inter-hemispheric 
frontal connections in ASD. Our volumetric finding is consistent 
with previous studies reporting reduced area of the entire corpus 
callosum, with greater magnitude of reduction in its anterior 
region (see the meta-analysis by Frazier and Hardan, 2009). 
Reduced FA in the corpus callosum is also consistent with many 
previous findings reported by others (Alexander et al., 2007; 
Noriuchi et al, 2010; Shukla et al, 2010). But, to the best of our 
knowledge, only one study reported a reduction in the number of 
inter-hemispheric frontal fibers using tractographic reconstruc- 
tions in patients with ASD (Thomas et al., 2011). Thomas et al. 
observed decreased numbers of streamlines specific to the body in 
high- functioning adults with ASD, which further correlated with 
ADI scores in the domain of restricted, repetitive and stereotyped 
behaviors. The fact that we replicate this correlation (though we 
focus on a more anterior, but overlapping, region) provides strong 
support for a role for the corpus callosum in repetitive behaviors, 
across ages and across IQ. 

REDUCED FA 

In addition to multimodal evidence for altered inter-hemispheric 
connectivity, we also observed seven clusters of decreased FA in 
the right hemisphere of patients with ASD compared to controls. 
The direction of our results was consistent with most previously 
published studies, which show decreased FA, although a few 
studies do report increased FA (reviewed in Travers et al., 2012). 



Surprisingly, in our small sample of children and adolescents with 
ASD, we found reduced FA only in the right hemisphere and did 
not detect changes to FA in the left hemisphere. Exclusively right- 
sided alterations to gyrification in the same sample of participants 
provide initial support for a relationship between white matter 
connectivity and cortical folding. However, it was not possible 
to detect whether the observed FA differences were related to 
differences in the degree of myelinisation or to differences in 
the orientation or number of white matter bundles, using voxel- 
wise measurements of FA. The spatial correspondence of altered 
gyrification and white matter microstructure in the same hemi- 
sphere led us to further examine the relationship between cortical 
folding and connectivity using more sophisticated tractographic 
measurements. 

CORRELATION BETWEEN GYRIFICATION AND CONNECTIVITY 

We did not observe a relationship between long-range connec- 
tivity and gyrification, as may have been expected from Van 
Essen's hypothesis that mechanical tension exerted on long con- 
nections shapes cortical folds (Van Essen, 1997). However, three 
out of the four clusters with decreased gyrification were mostly 
located in the frontal region, i.e., the region where an important 
decrease in inter-hemispheric connectivity was observed. The 
co-occurrence of decreased long-range connections in regions 
of altered gyrification points to a possible relationship between 
these two anatomical variables, but the mechanisms governing 
their association is likely to be more complex than what a linear 
regression can capture. 

We did, however, observe significant positive correlations 
between IGl and short-range connectivity in patients with ASD, 
but not in controls. This positive correlation means that higher 
IGl was observed in patients with higher intra-lobar (short-range) 
connectivity. According to Van Essen's theory, it may also be 
that short-range connections affect the creation of cortical folds 
during early brain development by reducing the distance between 
strongly interconnected regions from the two banks of one gyrus, 
thereby permitting compact wiring of the brain. Accordingly, the 
gyrification alterations observed in the present study may be a 
compensatory way of coping with altered connectivity in patients 
with ASD. 

LIMITATIONS AND CONCLUSION 

The main limitation of our study is its small sample size, restrict- 
ing our ability to identify age-related maturational changes or 
subtle brain-behavior relationships. We realize that, in a hetero- 
geneous disorder such as ASD, such small sample size may lead 
to observation of findings that may not be representative of the 
variability observed across the spectrum. However, despite the 
small sample size, we demonstrate the feasibility of multimodal 
studies in autism, bridging the gap between reports of altered 
cortical morphometry and findings of abnormal connectivity 
patterns. These preliminary results provide initial support for 
the idea that a higher degree of short-range connectivity alters 
the shape of the brain in patients with ASD during early neural 
development, and are an encouraging starting point for exploring 
this issue in larger samples of children, adolescents or adults with 
autism. 
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